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Gain and loss of chromosomal material is characteristic 
of bladder cancer, as wet! as malignant transformation in 
general. The consequences of these changes at both the 
transcription and translation levels is at present unknown 
partly because of technical limitations. Here we have at* 
tempted to address this question in pairs of non-Invasive 
and invasive human bladder tumors using a combination 
of technology that Included comparative genomic hybrid- 
ization, high density oligonucleotide array-based monitor- 
ing of transcript levels (5600 genes), and high resolution 



phenomenon at both the transcription and translation levels. 
High throughput array studies of the breast cancer cell line 
BT474 has suggested that there Is a correlation between 
DNA copy numbers and gene expression in highly amplified 
areas (2), and studies of individual genes In solid tumors 
have revealed a good correlation between gene dose and 
mRNA or protein levels in the case of c-erb~B2 f cyc//n d1 t 
emsl, and N-myc (3-5). However, a high cyclin D1 protein 
expression has been observed without simultaneous am- 



two-dimensional gel electrophoresisfthe results showed^ 11 * 103 * 10 " < 4 K and a low !eveI of c - m VC °°Py number in- 



that there is a gene dosage effect that in some cases 
superimposes on other regulatory mechanisms. This ef- 
fect depended (p < 0.015) on the magnitude of the com- 
parative genomic hybridization change. In general (18 of 
23 cases), chromosomal areas with more than 2-fold gain 
of DNA showed a corresponding increase in mRNA tran- 
scripts. Areas with loss of DNA, on the other hand, 
showed either reduced or unaltered transcript levels^; Be- 
cause most proteins resolved by two-dimensional gels 
are unknown it was only possible to compare mRNA and 
protein alterations in relatively few cases of well focused 
abundant proteins* 4vith few exceptions we found a good 
correlation (p < 0.005) between transcript alterations and 
protein levels. The implications, as well as limitations, 
of the approach are discussed. Molecular & Cellular 
Protoomlcs 1:37-45, 2002. 
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, Aneuploidy Is a common feature of most human cancers 
(1), but little is known about the genome-wide effect of this 
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crease was observed without concomitant c-myc protein 
overexpresslon (6). 

In human bladder tumors, karyotyping, fluorescent In situ 
hybridization, and comparative genomic hybridization (CGH) 1 
have revealed chromosomal aberrations that seem to be 
characteristic of certain stages of disease progression. In the 
case of non-Invasive pTa transitional cell carcinomas (FCCs), 
this includes loss of chromosome 9 or parts of it, as well as 
loss of Y in males. In minimally Invasive pT1 TCCs, the fol- 
lowing alterations have been reported: 2q-, 11p~, 1q+, 
11q13+, 17q+, and 20q+ (7-12). It has been suggested that 
these regions harbor tumor suppressor genes and onco- 
genes; however, the large chromosomal areas Involved often 
contain many genes, making meaningful predictions of the 
functional consequences of losses and gains yery difficult. 

In this investigation we^have combined gehome-wlde tech- 
nology for detecting genomic gains and losses (CGH) with 
gene expression profiling techniques (microarrays and pro- 
teomlcs) to determine the effect of gene copy number on 
transcript and protein levels In pairs of non-lnvaslve and In- 
vasive human bladder TCCs. 

EXPERIMENTAL PROCEDURES 

Material— Bladder tumor biopsies, were sampled after Informed 
consent was obtained and after removal of tissue for routine pathol- 
ogy examination. By llgjht microscopy tumors 335 and 532 were 
staged by an experienced pathologist as pTa (superficial papillary), 

1 The abbreviations used are: CGH, comparative genomic hybrid- 
ization; TCC, transitional cell carcinoma: UOH, ioss of heterozygosity; 
PA-FABP, psoriasis-associated fatty acid-binding protein; 2D, 
two-dimensional. 
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Fig. 1. DNAcopy number and mRNA expression level. Shown from left to right are chromosome (Cnr.), GGH profiles, gene location and 
expression level of specific genes, and overall expression level along the chromosome. A, expression of mRNA In Invasive tumor 733 as 
compared with the non-invasive counterpart tumor 335. B, expression of mRNA in invasive tumor 827 compared with the non-invasive 
counterpart tumor 532. The average fluorescent signal ratio between tumor DNA and normal DNA is shown along the length of the chromosome 
ifeft). The bold curve in the ratio profile represents a mean of four chromosomes and is surrounded by thin curves indicating one standard 
deviation. The central vertical fine broken) indicates a ratio value of 1 (no change), and the vertical fines next to It (dotted) indicate a ratio of 
0.5 (feft) and 2.0 (right). In chromosomes where the non-Invasive tumor 335 used for comparison showed alterations in DNA content, the ratio 
profile of that chromosome is shown to the right of the invasive tumor profile. The colored bars represents one gene each, Identified by the 
running numbers above the bars (the name of the gene can be seen at www.MDLDK/sdata,htmi). The bars Indicate the purported location of 
the gene, and the colors Indicate the expression level of the gene in the invasive tumor compared with the noninvasive counterpart: >2-fold 
Increase (WaeA), >2-fold decrease (o/ue), no significant change [orange). The bar to the far right, entitled Expression shows the resulting change 
In expression along the chromosome; the colors Indicate that at least half of the genes were^ up-regulated (Wac*), at least half of the genes 
dowivregulated (b/ue), or more than half of the genes are unchanged (orange). H a gene was absent intone of the samples and presenTTn 
another, it was regarded as more than a 2-fold change. A 2-fold level was chosen as this corresponded to one standard deviation In a double 
determination of -1800 genes. Centromeres and heterochromatic regions were excluded from data analysis. 



grade I and 11, respectively, tumors 733 and 827 were staged as pT1 
(invasive into submucosal 733 was staged as solid, and 827 was 
staged as papillary, both grade III. 

mRNA Preparation -Tissue biopsies, obtained fresh from surgery, 
were embedded Immediately In a sodium-guanidintum thiocyanate 
solution and stored at -80 °C. Total RNA was Isolated using the 
RNAzol B RNA isolation method (WAK-Chemle Medical GMBH). 
poly(A) + RNA was Isolated by an oOgo(oT) selection step (Ollgotex 
mRNA kit; Qiagen). 

cRNA Preparation-* fig of mRNA was used as starting material. 
The first and second strand cDNA synthesis was performed using the 
Superscripts choice system (Invitrogen) according to the manufac- 
turer's instructions but using an oRgo(cfT) primer containing a T7 RNA 
polymerase binding site. Labeled cRNA was prepared using the ME- 
GAscrip® in vitro transcription kft (Amblon). BiotMabeled CTP and 



UTP (Enzo) was used, together with unlabeled NTPs in the reaction. 
Following the in vitro transcription reaction, the unincorporated nu- 
cleotides were removed using RNeasy columns (Qiagen). 

Array Hybridization and Scanning— Array hybridization and scan* 
ning was modified from a previous method (13). 10 fig of cRNA was 
fragmented at 94 °C for 35 min In buffer containing 40 nut Tris 
acetate, pH 8.1 , 100 rrw KOAc, 30 nw MgOAc. Prior to hybridization, 
the fragmented cRNA In a 6x SSPE-T hybridization buffer (1 m NaCI, 
10 mM Tris, pH 7.6, 0.005% Triton), was heated to 95 °C for 5 min, 
subsequently cooled to 40 °C, and loaded onto the Affymetrix probe 
array cartridge. The probe array was then Incubated for 16 h at 40 °C 
at constant rotation (60 rpm). The probe array was exposed to 10 
washes in 6X SSPE-T at 25 °C followed by 4 washes in 0.5X SSPE-T 
at 50 *C. Trie biotinylated cRNA was. stained, with a streptavldln- 
phycoerythrin conjugate, 10 'jig/ml (Molecular Probes) In 6x SSPE-T 
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for 30 min at 25 °C followed by 1 0 washes in 6x SSPE-T at 25 °C. The 
probe arrays were scanned' at 560 nm using a confoca! laser scanning 
microscope {made for Affyrnetrix by Hewlett-Packard). The readings 
from the quantitative scanning were analyzed by Affyrnetrix gene 
expression analysis software. 

Microsatelfite Analysis — MJcrosatell He Analysis was performed as 
described previously (14). Microsatefiitea were selected by use of 
wwwj^br^Tm:nlh,g6v/genemap98 f and primer sequences were ob- 
tained from the genome data base at wwwgcfo.org. DNA was extracted 
from tumor and blood and amplified by PGR in a volume of 20 jxl for 35 
cycles. The ampficons were denatured and electrophoresed for 3 h in an 
ABI Prism 377. Data were collected tn the Gene Scan program for 
fragment analysis. Loss of heterozygosity was defined as less than 33% 
of one allele detected tn tumor ampiicons compared with bipod. 

Proteomfc Analysis— TCCs were minced into small pieces and 
homogenized in a small glass homogenizer In 0.5 ml of lysis solution. 
Samples were stored at -20 °C until use. The procedure for 2D gel 
electrophoresis has been described in detail elsewhere (15, 16), Gels 
were stained with sliver nitrate and/or Coornassie Brilliant Blue. Pro- 
teins were Identified by a combination of procedures that Included 
microsequenclng, mass spectrometry, hwp-drmenslonai gel Western 
immunoblotting, and comparison with the master two-dimensional gel 
image of human keratlnocyte proteins; see blobase.dk/cgl-bln/cells. 

CGH~Hybridlzatlon of differentially labeled tumor and normal DNA 
to normal metaphase chromosomes was performed as described 
previously (10). Ruoresceln-labeled tumor DNA (200 ng), Texas Red- 



labeled reference DNA (200 ng), and human Cot-1 DNA (20 /ig) were 
denatured at 37 °C for 5 mln and applied to denatured normal met- 
aphase slides. Hybridization was at 37 °C for 2 days. After washing, 
the slides were counterstained with 0.15 pg/ml 4,6-diamidino-2-prie- 
nyisndoie In an anti-fade solution. A second hybridization was per- 
formed for all tumor samples using fluoresceln-tabeted reference DNA 
and Texas Red-labeled tumor DNA (inverse labeling) to confirm the 
aberrations detected during the initial tiybrldizatton. Each CGH ex- 
periment also Included a normal control hybridization using fluores- 
cein- and Texas Red-labeled normal DNA. Digital image analysis was 
used to identify chromosomal regions with abnormal fluorescence 
ratios, Indicating; regions of DNA gains and losses. The average 
greenrred fluorescence intensity ratio profiles were calculated using 
four Images of each chromosome (eight chromosomes total) with 
normalization of the greercred fluorescence Intensity ratio for the 
entire metaphase and background correction. Chromosome identifi- 
cation was performed based on 4.6-diamldIno-2-phenylindole band- 
ing patterns. Only Images showing uniform high Intensity fluores- 
cence with minimal background staining were analyzed. All 
centromeres, p arms of acrocentric chromosomes, and heterochro- 
matic regions were excluded from the analysis. 

RESULTS 

Comparative Genomic Hybridization— Th* CGH analysis 
identified a number, of chromosomal gains and losses in the 
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Table I 

Correlation between alterations detected by COM and by expression monitoring 

Top, CGH used as Independent variable (If CGH alteration - what expression ratio was found); bottom, altered expression used as 
Independent variable (if expression alteration - what CGH deviation was found). 



CGH alterations 



Tumor 733 vs. 335 
Expression change clusters 



Concordance CGH alterations 



Tumor 827 vs. 532 
Expression change clusters 



Concordance 



13 Gain 
to I 



10 Up-regulation 

0 Down-regulation 

3 No change 

1 Up-regulation 

5 Down-regulation 

4 No change 



77% 



50% 



10 Gain 



12 Loss 



8 Up-regulation 
0 Down-regulation 

2 No change 

3 Up-regulation 

2 Down regulation 
7 No change 



80% 



17% 



Expression change clusters 



Tumor 733 vs. 335 
CGH alterations 



Concordance Expression change clusters 



Tumor 827 vs. 532 
CGH alterations 



Concordance 



16 Up-regulation 
21 Down-regulation 
15 No change 


11 Gain 

2 Loss 

3 No change 
1 Gain 

8 Loss 

12 No change 
3 Gain " 

3 Loss 

9 No change 


69% 
38% 
60% 


17 Up-regulation 
9 Down^regulation 
21 No change 


10 Gain 

5 Loss 

2 No change 
0 Gain 

3 Loss 

6 No change 
IGain 


59% 
33% 
81% 








3 Loss 

17 No change 





two Invasive tumors (stage pT1, TCCs 733 and 827), whereas 
the two non-Invasive papillomas (stage pTa, TCCs 335 and 
532) showed only 9p- f 9q22-q33-, and and 7+, 9q-, 
and Y— . respectively. Both invasive tumors showed changes 
(1q22-24+, 2q14.1-qter-, 3q12-q13.3-, 6q12-q22~, 
9q344, 11q12-q13+, 17+, and 20q11.2-q12+) that are typ- 
ical for their disease stage, as well as additional alterations, 
some of which are shown in Fig. 1. Areas with gains and 
losses deviated from the normal copy number to some extent, 
and the average numerical deviation from normal was 0,4-fold 
In the case of TCC 733 and 0.3-fold for TCC 827. The largest 
changes, amounting to at least a doubling of chromosomal 
content, were observed at 1q23 In TCC 733 (Fig, 1A) and 
20q12 in TCC 827 (Fig. 1B). 

mRNA Expression In Relation to DNA Copy Number— The 
mRN/yevels from the two lnvasiveJtumors w fTCCs 827 and 
733) were compared with the two non-invasive counterparts 
(TCCs 632 and 335). This was done in two separate experi- 
ments in which we compared TCCs 733 to 335 and 827 to 
532, respectively, using two different scaling settings for the 
arrays to rule out scaling as a confounding parameter. Ap- 
proximately 1,800 genes that yielded a signal on the arrays 
were searched In the Unigene and Genemap data bases for 
chromosomal location, and those with a known location 
(1096) were plotted as bars covering their purported locus. In 
that way It was possible to construct a graphic presentation of 
DNA copy number and relative mRNA levels along the indi- 
vidual chromosomes i)« 

For each mRNA a ratio was calculated between the level in 
the Invasive versus the non-Invasive counterpart Bars, which 
represent chromosomal location of a gene, were color-coded 
according to the expression ratio, and only differences larger 



than 2-fold were regarded as informative (Fig. 1). The density 
of genes along the chromosomes varied, and areas contain- 
ing only one gene were excluded from the calculations. The 
resolution of the CGH method is very tow, and some of the 
outlier data may be because of the fact that the boundaries of 
the chromosomal aberrations are not known at high resolution. 

Two sets of calculations were made from the data. For the 
first set we used CGH alterations as the independent variable 
and estimated the frequency of expression alterations in these 
chromosomal areas. In general, areas with a strong gain of 
chromosomal material contained a cluster of genes having 
increased mRNA expression. For example, both chromo- 
somes 1q21 -q25, 2p and 9q, showed a relative gain of more 
than 100% in DNA copy number that was accompanied by 
increased mRNA expression levels in the two tumor pairs (Fig. 
1). In most cases, chromosqrnai gains detected by CGH wer/j 
accompanied by an increased level of transcripts in both 
TCCs 733 (77%) and 827 (80%) (Table I, top). Chromosomal 
losses, on the other hand, were not accompanied by de- 
creased expression in several cases, and were often regis- 
tered as having unaltered RNA levels (Table I, top). The inabil^ 
ity to detect RNA expression changes In these cases was not 
because of fewer genes mapping to the lost regions (data not 
shown). 

In the second set of calculations we selected expression 
alterations above 2-fold as the Independent variable and es- 
timated the frequency of CGH alterations in these areas. As 
above, we found that Increased transcript expression corre- 
lated with gain of chromosomal material (TCC 733, 69% and 
TCC 827, 59%), whereas reduced expression was often de- 
tected In areas with unaltered CGH ratios (Table I, bottom). 
Furthermore, as a control we looked at areas with no after- 
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detected 
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Expression changes 
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Tumor 827 versus 532 Tumor 733 versus 335 

Ra 2. Correlation between maximum CGH aberration and the ability to detect expression change by oligonucleotide array 
monitoring. The aberration ts shown as a numerical -fold change In ratio between invasive tumors 827 (A) and 733 (4 ) arid their non-Invasive 
counterparts 532 and 335. The expression change was taken from the Expression line to the right In Fig. 1, which depicts the resulting 
expression change for a given chromosomal region. At least half of the mRNAs from a given region have to be either up- or down-regulated 
to be scored as an expression change. All chromosomal arms In which the CGH ratio plus or minus one standard deviation was outside the 
ratio value of one were included. 



ation in expression. No alteration was detected by CGH In. 
most of these areas (TCC 733, 60% and TCC 827, 81 %; see 
fable I, bottom). Because the ability to observe reduced or 
Increased mRNA expression clustering to a certain chromo- 
somal area clearly reflected the extent of copy number 
changes, we plotted the maximum CGH aberrations In the 
regions showing CGH changes against the ability to detect a 
change In mRNA expression as monitored by the oligonucleo- 
tide arrays (Fig. 2)(B>r both tumors TCC 733 (p < 0.015) and 
TCC 827 (p < 0.00003) a highly significant correlation was 
observed between the level of CGH ratio change (reflecting 
the DNA copy number) and alterations detected by the array 
based technology (Rg. 2^ Similar data were obtained when 
areas with altered expression were used as Independent vari- 
ables. These areas correlated best with CGH when the CGH 
ratio deviated 1.6- to 2.0-fold (Table I, bottom) but mostly did 
not at lower CGH deviations. These data probably reflect that 
loss of an allele may only lead to a 50% reduction In expres- 
sion level, which is at the cut-off point for detection of expres- 
sion alterations. Gain of chromosomal material can occur to a 
much larger extent, 

Mlcrvsatellfte-based Detection of Minor Areas of Loss- 
es—In TCC 733, several chromosomal areas exhibiting DNA 
amplification were preceded or followed by areas with a nor- 
mal CGH but reduced mRNA expression {see Fig. 1 , TCC 733 
chromosome 1q32, 2p21, and 7q21 and q32, 9q34, and 
10q22). To determine whether these results were because of 
undetected loss of chromosomal material In these regions or 



because of other non-structural mechanisms regulating tran- 
scription, we examined two mlcrosatellltes positioned at chro- 
mosome 1q25-32 and two at chromosome 2p22. Loss of 
heterozygosity (LOH) was found at both 1q25 and at 2p22 
indicating that minor deleted areas were not detected with the 
resolution of CGH (Rg. 3). Additionally, chromosome 2p in 
TCC 733 showed a CGH pattern of gain/no change/gain of 
DNA that correlated with transcript increase/decrease/in- 
crease. Thus, for the areas showing increased expression 
there was a correlation with the DNA copy number alterations 
(Fig. 1A). As indicated above, the mRNA decrease observed in 
the middle of the chromosomal gain was because of LOH, 
implying tri&t one of the mechanism^ for mRNA down-regu- 
lation may be regions that have undergone smaller losses of 
chromosomal material. However, this cannot be detected with 
the resolution of the CGH method. 

In both TCC 733 and TCC 827, the telomeric end of chro- 
mosome 11p showed a normal ratio In the CGH analysis; 
however, clusters of five and three genes, respectively, lost 
their expression. Two mlcrosatellltes (D11S1760, D11S922) 
positioned close to MUC2, IGF2, and cathepsln D Indicated 
LOH as the most likely mechanism behind the loss of expres- 
sion (data not shown). 

A reduced expression of mRNA observed in TCC 733 at 
chromosomes 3q24> 11p11, 12p12.2, 12q21.1, and 16q24 
and in TCC 827 at chromosome 11p15.5, 12p11, 15q11.2, 
and 18q12 was also examined for chromosomal losses using 
mlcrosatellltes positioned as close as possible to the gene loci 
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Fig. 3. Microsateilrte analysis of loss of heterozygosity. Tumor 
733 showing loss of heterozygosity at chromosome 1q25, detected 
(a) by D1S215 close to Hu class I histocompatibility antigen (gene 
number 38 In Rg. 1), (b) by D1S2735 close to cathepsin E (gene 
number 41 In Fig. U and (c) at chromosome 2p23 byja2§^2$4. close 
to general 0-spectrin (gene number 11 on Rg. 1) and of (d) tumor 827 
showing loss of heterozygosity at chromosome 18q12 by S18S1118 
close to mitochondrial 3-oxoacyl-coenzyrne A thiolase (gene number 
12 In Rg. 1). The upper curves show the electropherogram obtained 
from normal DNA from leukocytes (A/), and the lower curves show the 
electropherogram from tumor DNA (7). In all cases one allele is 
partially lost In the tumor ampllcon. 

showing reduced mRNA transcripts. Only the mlcrosatellite 
positioned at 18q12 showed LOH (Fig. 3), suggesting that 
transcriptional down-regulation of genes In the other regions 
may be controlled by other mechanisms. 

Relation between Changes In mRNA and Protein Levels- 
2D-PAGE analysis, in combination with Coomassle Brilliant 
Blue and/or silver staining, was carried out on all four tumors 
using fresh biopsy material, 40 well resolved abundant known 
proteins migrating In areas away from the edges of the pH 



A B 

Reduced Unaltered Increased 2*0~get 

protein protein protein 

Ra 4. Correlation between protein levels as Judged by 2D* 
PAGE and transcript ratio. For comparison proteins were divided In 
three groups, unaltered In level or up- or down-regulated horizontal 
axis). The mRNA ratio as determined by oligonucleotide arrays was 
plotted for each gene [vertical axis). A, mRNAs that were scored as 
present In both tumors used for the ratio calculation; A, mRNAs that 
were scored as absent in the invasive tumors {along horizontal axis} or 
as absent In non-Invasive reference {top of figure). Two different 
scallngs were used to exclude scaling as a confounder, TCCs 827 
and 532 (AA) were scaled with background suppression, and TCCs 
733 and 335 (#0) were scaled without suppression. Both compari- 
sons showed highly significant (p < 0.005) differences in mRNA ratios 
between the groups. Proteins shown were as follows: Qroup A (from 
fe/r), phosphoglucomutase 1, glutathione transferase class jt number 
4, fatty acid-binding protein homologue, cytokeratln 15, and cyto- 
keratln 13; B (from left), fatty acld-blndlng protein homologue, 28-kDa 
heat shock protein, cytokeratln 13, and calcyclin; C<from terr), a-eno- 
lase, hnRNP 81. 28-kDa heat shock protein, 14-3-3-e, and 
pre-mRNA splicing factor; 0, mesothelial keratin K7 (type If); E (from 
fop), glutathione S-transferase-Tr and mesothelial keratin K7 (type. 10; 
F (from top and left), adenytyl cyclase-associated protein, E-cadherin, 
keratin 19, calgizzarin, pnosphoglycerate mutase, . annexln IV, cy- 
toskeletal -y^actln, hnRNP AT, Integral membrane protein calnetxln 
(IP90), hnRNP H, brain-type clathrin light chain-a, hnRNP F, 70-kDa 
heat shock protein, heterogeneous nuclear ribonudeoprotein A/B, 
translationaily controlled tumor protein, liver gIyceraldehyde-3-phos- 
phate dehydrogenase, keratin 8, aldehyde reductase, and Na,K- 
ATPase 0-1 subunit; G, (from top and fe$, TCP20, calgizzarin. 70- 
kDa heat shock protein, calnexln, hnRNP K, cytokeratirt 15, ATP 
synthase, keratin 1 9, trlosephosphate feomerase, hnRNP F, liver glyc- 
eraldehyde-3 -phosphatase dehydrogenase, glutathione S-transfer- 
ase-Tr, and keratin 8; H (from fefl), plasma gelsdlin, autoantigen cal- 
retlculin, thioredoxin, and NAD+ -dependent 15 hydroxyprostaglandin 
dehydrogenase; I (from top), prolyl 4-hydroxylase 0-subunlt, cyto- 
keratin 20, cytokeratln 17, prohibition, and fructose triphos- 
phatase; J annexln II; K, annexln IV; L (from top and left), 90-kDa heat 
shock protein, prolyl 4-hydroxyiase 0-subunit, a-enolase, GRP 78, 
cyclophilin, and cofilln. 

gradient, and having a known chromosomal location, were 
selected for analysis in the TCC pair 827/532. Proteins were 
identified by a combination of methods (see "Experimental 
Procedures 0 ). In general there was a highly significant corre- 
lation (p < 0,005) between mRNA and protein alterations (Fig. 
4). Only one gene showed disagreement between transcript 
alteration and protein alteration. Except for a group of cyto- 
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keratins encoded by genes on chromosome 17 (Fig. 5) the 
analyzed proteins did not belong to a particular f amity. 26 well 
focused proteins whose genes had a know chromosomal 
location were detected in TCCs 733 and 335, and of these 19 
correlated (p < 0,005) with the mRNA changes detected using 
the arrays (Fig. 4). For example, PA-FABP was highly ex- 
pressed in the non-invasive TCC 335 but lost in the Invasive 
counterpart (TCC 733; see Fig. 5). The smaller number of 
proteins detected in both 733 and 335 was because of the 
smaller size of the biopsies that were available. 

11 chromosomal regions where CGH showed aberrations 
that corresponded to the changes in transcript levels also 
showed corresponding changes in the protein level (Table II). 
These regions Included genes that encode proteins that are 
found to be frequently altered in bladder cancer! namely 
cytokeratins 17 and 20, annexins II and IV, and the fatty 
acid-binding proteins PA-FABP and FBP1. Four of these pro- 
teins were encoded by genes In chromosome 17q f a fre- 
quently amplified chromosomal area In invasive bladder 
cancers. 

DISCUSSION 

Most human cancers have abnormal DNA content, having 
lost some chromosomal parts and gained others. The present 
study provides some evidence as to the effect of these gains 
and losses on gene expression in two pairs of nqn-lnyaslye 
and Invasive TCCs using high throughput expression arrays 
and proteomics, in combination with CGH. In general, the 
results showed that there is a clear individual regulation of the 
mRNA expression of single genes, which in some cases was 
superimposed by a DNA copy number effect Inmost cases, 
genes located In chromosomal areas with gains often exhib- 
ited increased mRNA expression, whereas areas showing 
losses showed either no change or a reduced mRNA expres- 
sion. The latter might be because of the fact that losses most 
often are restricted to loss of one allele, and the cut-off point 
for detection of expression alterations was a 2-foW change, 
thus being at the border of detection. In several cases, how- 

.Table II 

Proteins whose expression level correlates with both mRNA and gene dose changes 



Protein Chromosomal location Tumor TCC CGH alteration Transcript alteration* Protein alteration 
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17q12-q21 


827 


Gain 


3.8-Fold up 


Increase 


Cytokeratin 20 


17q21.1 


827 


Gain 


5.6-Fold up 


Increase 


(PA-)FABP 


8q21.2 


827 


Loss 


10-Fold down 


Decrease 


FBP1 


9q22 


827 


Gain 


2.3-Fold up 


Increase 


Plasma gelsolin 


9q31 


827 


Gain 


Abs to Pres 


Increase 


Heat shock protein 28 


15q12-q13 


827 


Loss 


2,5-Fold up 


Decrease 


Prohibltin 


17q21 


827/733 


Gain 


3.7V2.5-Fold up* 


Increase 


Prolyl-4-hydroxyl 


17q25 


827/733 


Gain 


5.7-/1 .6-Fold up 


Increase 


hnRNPBt 


7p15 


827 


Loss 


2.5-Fold down 


Decrease 



a Abs, absent; Pres, present. 

* In cases where the corresponding alterations were found In both TCCs 827 and 733 these are shown as 827/733. 




Fig. 5. Comparison of protein and transcript levels in invasive 
and non-Invasive TCCs, The upper part of the figure shows a 2D gel 
ifeft) and the oligonucleotide array Iftghfj of TCC 532. The red rectan- 
gles on the upper get highlight the areas that are compared below. 
Identical areas of 2D gels of TCCs 532 and 827 are shown below. 
Clearly, cytokeratins 13 and 15 are strongly down-regulated In TCC 
827 (reef annotation). The tile on the array containing probes for 
cytokeratin 15 Is enlarged below the array (red arrow) from TCC 532 
and Is compared with TCC 827. The upper row of squares in each tile 
corresponds to perfect match probes; the lower row corresponds to 
mismatch probes containing a mutation (used for correction for un- 
specific binding). Absence of signal is depicted as black, and the 
higher the signal the lighter the color. A high transcript level was 
detected In TCC 532 (6151 units) whereas a much lower level was 
detected in TCC 827 (absence of signals). For cytokeratin 13, a high 
transcript lever was also present in TCC 632 (15659 units), and a 
much lower level was present in TCC 827 (623 units). The 2D gels at 
the bottom of the figure QefQ show levels of PA-FABP and adipocyte- 
FABP in TCCs 335 and 733 (invasive), respectively. Both proteins are 
down-regulated in the invasive tumor. To the right we show the array 
tiles for the PA-FABP transcript A medium transcript level was de- 
tected In the case of TCC 335 (1 277 units) whereas very low levels 
were detected In TCC 733 (1 66 units). IEF, Isoelectric focusing. 



Molecular & Cellular Proteomlcs 1.1 43 



Qene Copy Numbers, Transcripts, and Prtrtein Levels 



ever, an increase or decrease in DNA copy number was 
associated with ere novo occurrence or complete loss of tran- 
script, respectively. Some of these transcripts could not be 
. detected In the non-invasive tumor but were present at rela- 
tively high levels In areas with DNA amplifications in the inva- 
sive tumors (e.g. in TCC 733 transcript from cellular ligand of 
annexin II gene (chromosome 1q21) from absent to 2670 
arbitrary Units; in TCC 627 transcript from small proline-rich 
protein 1 gene (chromosome 1q12~q21.l) from absent to 
1326 arbitrary units). It may be anticipated from these data 
that significant clustering of genes with an increased expres- 
sion to a certain chromosomal area indicates an increased 
likelihood of gain of chromosomal material in this area. 

Considering the many possible regulatory mechanises act- 
ing at the level of transcription, it seems striking that the gene 
dose effects were so clearly detectable in gained areas. One 
hypothetical explanation may lie In the loss of controlled 
methylation In tumor cells (17-19). Thus, it may be possible 
that Iri chromosomes with increased DNA copy numbers two 
. or more alleles could be demethyiated simultaneously leading 
to a higher transcription level, whereas In chromosomes with 
losses the remaining allele could be partly methylated, turning 
off the process (20, 21). A recent report has documented a 
ploidy regulation of gene expression in yeast, but in this case ail 
the genes were present In the same ratio (22), a situation that is 
not analogous to that of cancer ceils, which show nrarked 
chromosomal aberrations, as well as gene dosage effects^ 

Several CGH studies of bladder cancer have shown that 
some chromosomal aberrations are common at certain 
stages of disease progression, often occurring in more than 1 
of 3 tumors. In pTa tumors, these include 9p- t 9q-, 1 q+, Y- 
(2, 6), and in-pT1 tumors, 2q-.11p-.11q-, 1q+, 5p+, 8q+, 
17q+, and 20q+ (2-4, 6, 7). The pTa tumors studied here 
showed similar aberrations such as 9p- and 9q22-q33- and 
9q- and Y-, respectively. Likewise, the two minimal invasive 
pT1 tumors showed aberrations that are commonly seen at 
that stage, and TCC 827 had a remarkable resemblance to the 
commonly seen pattern ofOosses and gains, such as 1q2g^24 
amplification (seen in both tumors), 1 1q14-q22 loss, the latter 
often linked to 1 7 q+ (both tumors), and 1q+ and 9p-, often 
linked to 20q+ and 11 q13+ (both tumors) (7-9). These ob- 
servations indicate that the pairs of tumors used In this study 
exhibit chromosomal changes observed in many tumors, and 
therefore the findings could be of general importance for 
bladder cancer. 

Considering that the mapping resolution of CGH Is of about 
20 megabases it is only possible to get a crude picture of 
chromosomal instability using this technique. Occasionally, 
we observed reduced transcript levels close to or inside re- 
gions with Increased copy numbers. Analysis of these regions 
by positioning heterozygous microsateliites as close as pos- 
sible to the locus showing reduced gene expression revealed 
loss of heterozygosity in several cases. It seems likely that 
multiple and different events occur along each chromosomal 



arm and that the use of cDNA microarrays for analysis of DMA 
copy number changes will reach a resolution that can resolve 
these changes, as has recently been proposed (2). The outlier 
data were not more frequent at the boundaries of the CGH 
aberrations. At present we do not know the mechanism be- 
hind chromosomal aneuploldy and cannot predict whether 
chromosomal gains will be transcribed to a larger extent than 
the two native alleles. A mechanism as genetic imprinting has 
an Impact on the expression level in normal ceils and is often 
reduced in tumors. However, the relation between imprinting 
and gain of chromosomal material is not known. 

We regard it as a strength of this Investigation that we were 
able to compare invasive tumors to benign tumors rather than 
to normal urothelium, as the tumors studied were biologically 
very dose, and probably may represent successive steps in 
the progression of bladder cancer. Despite the limited amount 
of fresh tissue available it was possible to apply three different 
state of the art methods. The observed correlation between 
DNA copy number and mRNA expression is remarkable when 
one considers that different pieces of the tumor biopsies were 
used for the different sets of experiments. This indicate that 
bladder tumors are relatively homogenous, a notion recently 
supported by CGH and LDH data that showed a remarkable 
similarity even between tumors and distant metastasis (10, 23). 

In the few cases analyzed, mRNA. and protein levels 
showed a striking correspondence although in some cases 
we found discrepancies that may be attributed to translatlonal 
regulation, post-translatlonal processing, protein degrada- 
tion, or a combination of these. Some transcripts belong to 
under-translated mRNA pools, which are associated with few 
translationally inactive ribosomes; these pools, however, 
seem to be rare (24). Protein degradation, for example, may 
be very Important in the case of polypeptides with a short 
half-life (e.g. signaling proteins). A poor correlation between 
mRNA and protein levels was found in liver cells as deter- 
mined by arrays and 2D^PAGE (25), and a moderate correla- 
tion was recently reported by Ideker et al. (26) In yeast 
. (Interestingly, our study revealed a much better correlation 
between gained chromosomal areas and increased mRNA 
levels than between loss of chromosomal areas and reduced 
mRNA levels. In general, the level of CGH change determined 
the ability to detect a change in transcript) One possible 
explanation could be that by losing one allele the change in 
mRNA level is not so dramatic as compared with gain of 
material, which can be rather unlimited and may lead to a 
severalfold Increase in gene copy number resulting In a much 
higher impact on transcript level, the latter would be much 
easier to detect on the expression arrays as the cut-off point 
was placed at a 2-fold level so as not to be biased by holse on 
the array. Construction of arrays with a better signal to noise 
ratio may In the future allow detection of lesser than 2-fold 
alterations in transcript levels, a feature that may facilitate the 
analysis of the effect of loss of chromosomal areas on tran- 
script levels. 
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: In eleven cases we found a significant correlation between 
. DMA copy number, mRNA expression; and protein level Four 
of these proteins were encoded by genes located at a fre- 
quent!/ amplified area in chromosome 17q.. Whether DNA 
copy number is one of the mechanisms behind alteration of 
these eleven proteins is at present unknown and will have to 
be proved by other methods using a larger number of sam- 
ples. One factor making such studies complicated Is the large 
extent of protein modification that occurs after translation, 
requiring immunoidentffication and/or mass spectrometry to 
correctly Identify the proteins In the gels. 

In conclusion, the results presented In this study exemplify 
the large body of knowledge that may be possible to gather in 
the future by combining state of the art techniques that follow 
the pathway from DNA to protein (26). Here, we used a tradi- 
tional chromosomal CGH method, but in the future high reso- 
lution CGH based on microarrays with many thousand radiation 
hybrid-rriapped genes will increase the resolution and informa- 
tion derived" from these types of experiments (2). Combined with 
expression arrays analyzing transcripts derived from genes with 
known locations, and 2D gel analysis to obtain Information at 
the post-translational level, a clearer and more developed un- 
derstanding of the tumor genome will be forthcoming. 
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